      SUBROUTINE SETUP
      RETURN
      END
C   HS42
      BLOCK DATA
      INCLUDE 'BLOCO.INC'
      INTEGER I,J
      DATA NAME/'HS42'/
      DATA (X(I),I=1,4)/4*1.0D0/
      DATA N/4/ , NH/2/ , NG/ 0/
      DATA DEL0/1.00D0/ ,TAU0/1.00D0/ ,TAU/.1D0/
      DATA ((GUNIT(I,J),I=1,3),J=0,2)/-1,0,0,-1,0,0,-1,0,0/
      END
      SUBROUTINE EF(X,FX)
      INCLUDE 'FUCO.INC'
      DOUBLE PRECISION X(NX),FX
      ICF=ICF+1
      FX=(X(1)-1.D0)**2 + (X(2)-2.D0)**2+(X(3)-3.D0)**2+(X(4)-4.D0)**2
      RETURN
      END
      SUBROUTINE EGRADF(X,GRADF)
      INCLUDE 'FUCO.INC'
      DOUBLE PRECISION X(NX),GRADF(NX)
      ICGF=ICGF+1
      GRADF(1)=2.D0*(X(1)-1.D0)
      GRADF(2)=2.D0*(X(2)-2.D0)
      GRADF(3)=2.D0*(X(3)-3.D0)
      GRADF(4)=2.D0*(X(4)-4.D0)
      RETURN
      END
      SUBROUTINE EH(I,X,HXI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),HXI
      CRES(I)=CRES(I)+1
      IF(I .EQ. 2)      GOTO 100
      HXI=X(1)- 2.D0
      RETURN
  100 HXI=X(3)**2+X(4)**2-2.D0
      RETURN
      END
      SUBROUTINE EGRADH(I,X,GRADHI)
      INCLUDE 'FUCO.INC'
      INTEGER I,K
      DOUBLE PRECISION X(NX),GRADHI(NX)
      CGRES(I)=CGRES(I)+1
      DO      10      K=1,4
      GRADHI(K)=0.D0
  10  CONTINUE
      IF(I .EQ. 2)      GOTO 100
      GRADHI(1)= 1.D0
      RETURN
  100 CONTINUE
      GRADHI(3)= 2.D0*X(3)
      GRADHI(4)= 2.D0*X(4)
      RETURN
      END
      SUBROUTINE EG(I,X,GXI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),GXI
      RETURN
      END
      SUBROUTINE EGRADG(I,X,GRADGI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX) ,GRADGI(NX)
      RETURN
      END

